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We calculate charmed meson spectra and decay constants in lattice QCD employing one-loop 0{a) improved 
heavy quark action and axial- vector currents. In quenched simulations at a ~ 0.1 fm with the plaquette gauge 
action as well as a renormalization-group improved one, it is shown that the deviation from the continuum 
dispersion relation and the violation of space-time symmetry for the pseudoscalar meson decay constants are 
substantially reduced, once the 0{a) improvement is applied. Preliminary results with two flavors of dynamical 
quarks are also presented. 



1. Introduction 

With current computational resources it is 
difficult to achieve a precise determination of 
weak matrix elements associated with B and 
D mesons because of large mga corrections for 
b- and c-quarks. Recently a relativistic heavy 
quark(RHQ) action was proposed to control 
these mga corrections by applying the on-shell 
0{a) improvement program to massive quarks^, 
and the mass-dependent improvement coefficients 
have been perturbatively determined up to one- 
loop level|2]. The improvement of the heavy- 
heavy and heavy-light axial vector currents are 
also carried out up to one-loop level in a mass 
dependent wayjjj. In this report we perform a 
numerical study of the RHQ action in quenched 
QCD by investigating the effective speed of light 
Coff from the meson energy and the space-time 
asymmetry R for the pseudoscalar decay con- 
stant. The deviation from Ccs = 1 and R = I 
measures typical systematic errors due to finite 
mqa corrections left in the formulation, which 
should be incorporated in the results of charmed 



meson spectra and decay constants. Preliminary 
results with two flavors of dynamical quarks are 
also presented. 

2. Simulation details 

We employ the clover(CL) action with non- 
perturbative Cg*^ |4I5| for light quarks and the 
RHQ action for heavy quarks combined with the 
plaquette and the Iwasaki gauge actions. It is 
worth mentioning that we replace the massless 
part of the one-loop improvement coefficients Cb 
and C£; in the RHQ action with that of the non- 
pcrturbative one as 
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where the superscript PT represents "perturba- 
tive" value. With this replacement 0{a) errors 
are completely removed at mg = 0. We calculate 
charmed meson spectra using 300 configurations 
on a 24"^ x 48 lattice at a ^ O.lfm in quenched 
QCD. Three values of the lihgt quark mass cor- 
responding to Mps/Mv ^ 0.56 — 0.77 and four 
values of the heavy quark mass around the charm 
quark mass are adopted. For comparison we per- 
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Figure 1. Effective speed of light for the PS me- 
son(top) and space-time asymmetry for the PS 
meson decay constant (bottom) as a frmction of 
the meson mass in quenched QCD. 
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Figure 2. Charmonium S-state hyperfine spht- 
ting(top) and Dg meson decay constant (bottom) 
in quenched QCD as a function of a set by Som- 
mer scale vq. 



form parallel calculations with the CL action for 
both heavy and light quarks at the same gauge 
coupling and lattice size. 

3. Numerical study of the RHQ action 

Let us examine effectiveness of the 0{a) im- 
provement for the RHQ action in a nonperturba- 
tive way. We extract the effective speed of light 
Ceff by fitting the pesudoscalar(PS) meson energy 
aE as a function of spatial momentum aps'. 

{aEf = (aMpoic)' + cl^^iapsf, 

where Mpoie denotes the pole mass. Ccs is ex- 
pected to become unity in the continuum limit. 
Numerical results are plotted in Fig. H^top) as a 
function of the meson mass for heavy-heavy(H-H) 
and heavy- light (H-L) systems. Around charmed 
meson mass the deviation from Coff = 1 is equal to 
or larger than about 10% for the CL heavy quark 
action, while less than 5% for the RHQ action. 

We also measure the space-time asymmetry in 
terms of the PS meson decay matrix elements: 

{Q\A^\PS) E_ 
^ - {i)\A^\PS) ■ 

which should become unity in the continuum 
limit. Results are plotted in FigQIbottom) as 



a function of the meson mass. The asymmetry 
reaches 10 ^ 20% aroimd the charmed meson 
mass for the CL action, while it is less than 8% 
for the RHQ action. These observations allow us 
to conclude that the improved action works as 
designed in reducing the mqa corrections. 

Figure [5Ja) shows the charmonium S-state hy- 
perfine splitting as a function of lattice spacing, 
where we practice two analysis: one is based on 
the pole mass of the charmed meson and another 
on the kinetic mass: 

aA/kin = aAfpoio/c^ff. 

We find 15% discrepancy, which should be in- 
cluded as a systematic error. Compared to re- 
sults with anisotropic lattice QCD[H|, our results 
are slightly larger at a ~ O.lfm. Both are sub- 
stantially lower than the experimental value. 

In Fig. Elb) we plot the Dg meson decay con- 
stant which are obtained from the time and spa- 
tial components of the axial vector current. The 
difference between two definitions is less than 3% 
for both gauge actions. For comparison we also 
plot the results obtained by the CL heavy quark 
action with the plaquette gauge action 0. Our 
results at a ~ O.lfm are much closer to the value 
of the CL heavy quark action in the continuum 
limit than at a similar lattice spacing. This would 
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Figure 3. Effective speed of light for tlie PS me- 
son(top) and space-time asymmetry for the PS 
meson decay constant (bottom) as a function of 
the meson mass in lattice unit(left) and physical 
unit(right) for Nf^2 QCD. 



suggest that the RHQ action indeed improves the 
scaling behavior of the decay constant. 

4. Preliminary results of N f — 2 QCD 

Encouraging results in quenched QCD urge us 
to apply the calculation to the two-flavor QCD 
gauge configurations |HI. We employ the mean- 
field improved CL light quark action and the 
RHQ action for the valence quraks. Up to now 
two set of lattices have been analyzed: 400/200 
configurations for 16^ x 32/24^ x 48 at /3 = 
1.95/2.1 with four sea quark masses. Lattice 
spacings are roughly O.lSfm and O.llfm, respec- 
tively. We make a similar choice for the set of 
valence heavy and light quark masses as in the 
quenched case. 

In Fig. 131 we show the effective speed of light 
and the space-time asymmetry for the PS me- 
son decay constant. As a decreases, the devi- 
ation from Coff = 1 and i? = 1 is reduced at 
fixed mqa and mg. These behavior are consis- 
tent with the expected leading scaling violation of 
the RHS action given by 0(/2(r7iQa)(aAQCD)^) or 
0{fi{mQa)g'^ah.QCT)), where /i and /2 are some 
unkown functions. A little surprising is that, as 
seen from the left panels in Fig. 13 /1.2 look almost 
linear in terms of mqa even beyond mqa — 1. 

Figure 0] illustrates the a dependence of the Dg 
meson decay constant. While the scaling behav- 
ior from with the RHQ action is only a little 
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Figure 4. Dg meson decay constant for Nf = 2 
QCD as a function of a. The scale is set by tq for 
our results and mp for FNAL method. 



milder over the previous result obtained by the 
FNAL interpretation with the CL heavy quark 
action|3, the one from Ak is much improved. 
This result suggests that we will be able to make a 
reliable continuum extrapolation, employing a si- 
multaneous fit for /d, obtained from Ai^ and Ak, 
once we have data on CP-PACS configurations at 
/3 = 2.2 or 1.8. 
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